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O 1 

■ We consider a novel mechanism for J/^ production in nuclear collisions arising due to the high 

| density of gluons. The resulting J/'fy production cross section is evaluated as a function of rapidity 

and centrality. We compute the nuclear modification factor and show that the rapidity distribution 
of the produced J/^'s is significantly more narrow in AA collisions due to the gluon saturation 
effects. Our results indicate that gluon saturation in the colliding nuclei is a significant source of 
J/\P suppression and can explain the experimentally observed rapidity and centrality dependencies 
of the effect. 



=5 



X 



I. INTRODUCTION 



The mechanism of J/ production in high energy nuclear collisions can be expected to differ from that in hadron- 
hadron collisions. Consider first the J production in hadron-hadron collisions. The leading contribution is given 
by the two-gluon fusion (i) G + G — ► J/ty + soft gluon, see Fig. Q1A. This process is of the order 0{a b s ). The 
three-gluon fusion (ii) G + G + G — > J/^, see Fig. [TJB, is parametrically suppressed as it is proportional to 0{oP s ). 
However, in hadron-nucleus collisions an additional gluon can be attached to the nucleus. This brings in an additional 
■s^j- ■ enhancement by a factor ~ A 1 / 3 . If the collision energy is high enough, the coherence length becomes much larger than 
| the size of the interaction region. In this case all A nucleons interact coherently as a quasi-classical field [H, H, H, IH . 
In the quasi-classical approximation a 2 A 1//3 ~ 1. Therefore, the three-gluon fusion is actually enhanced by l/a s as 
compared to the two-gluon fusion process. Similar conclusion holds for heavy-ion collisions (we do not consider any 
qq ■ final state processes leading to a possible formation of the Quark-Gluon Plasma). 

In this letter we calculate the rapidity and centrality dependence of J/^f production in A A collisions taking into 
account the gluon saturation effects 0, Q. The case of J/^f production in dAu collisions was considered by two of us 
some time ago 0] where exactly the same mechanism was discussed, though the nuclear geometry was oversimplified. 
We note a reasonable agreement of this earlier approach with the dA data; a more definite conclusion can be reached 
once the higher statistics data becomes available in the near future. Our goal is to understand to what extent the cold 
nuclear matter effects are responsible for the J/\& suppression in AA collisions. In our calculation we use the dipole 
model [§[ and take into account realistic nuclear geometry. In Ref. 0] detailed arguments were given which justify the 
application of the dipole model to the calculation of J/ty production at RHIC. It was argued that (i) the coherence 
length for the production of the cc pair is sufficiently larger than the longitudinal extent of the interaction region. This 
means that the development of the light-cone "wave function" happens a long time before the collision, (ii) Formation 
of J is characterized by a time scale on the order of the inverse binding energy. This time is certainly much larger 
than the cc production coherence length (by a factor of ~ l/a 2 .) implying that the formation process takes place far 
away from the nucleus. Therefore, in the following, we concentrate on the dynamics of cc pair interactions with the 
nucleus. 



II. NEW MECHANISM OF J/* PRODUCTION 

A particularly helpful insight into the nature of the new production mechanism, depicted in Fig.[T]-B is obtained if we 
note that the three-gluon contribution is suppressed as compared to the two-gluon fusion mechanism (Fig. [ljA) by an 
additional factor r 2 , where r is cc transverse separation such that (2m c ) _1 < r < (2to c O! s ) . This factor arises since 
we need to have three gluons in the area of the order of r 2 . In other words, it means that this reaction is originated 
from the next-twist contribution. However, in the hadron - nucleus interactions the next-twist contribution appears 
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FIG. 1: The process of inclusive J/^ production in hadron-hadron (Fig. QJA) and in hadron-nucleus collisions (Fig. QJB). 

always in the dimensionless combination r 2 Q 2 s with the saturation scale Q s . The saturation scale is proportional to 
A 1 / 3 which compensates for the smallness of r. The dominance of the higher twist process is the main idea of Q and 
we develop it in this paper in the case of heavy-ion collisions. 

Fig. [1}B represents the contribution of the order (a 2 A 1 / 3 ) 2 . In general, there must be an odd number of gluons 
connected to the charm fermion line because the quantum numbers of J/\I / and of gluon are 1 . Therefore, each 
inelastic interaction of the cc pair must involve two nucleons and hence is of the order (a;: A 1 / 3 ) 2 ™, where n — 
1,2,..., A/2 is the number of nucleon pairs. To take this into account we write the cross section as the sum over all 
inelastic processes (labeled by the index n). This sum involves only even number of interactions. For a heavy nucleus 
A> 1 we have 

da in (pA) 2 
dY( pb = C F x ^ G \ x ^ m c) (!) 

f2R A r r 

pa m (x 2 ,r,r')dz / d 2 r V G (h, r, z = 1/2) V v (r) ® / d 2 r'V G {h, r' , z = 1/2) V v (r') 

oo r 2R A „2R A r 2R A \ 

(<r(i2 ,r 2 ) + <r(x 2 ,r' 2 )) P 2R A £ / dg l dz2 . . . / dz 2n+1 p 2n ^ a 2 ^ 1 (x 2 , r, r>)\ 

n=0"' z ° J z 2n ) 

where ^l G ® is the projection of the J/^/ light-cone "wave-function" onto the virtual gluon one 0,0], z%s are the 
nucleon longitudinal coordinates in the nucleus and 

Vin(x2,r,r') = a{x 2 ,r 2 ) + a(x 2l r' 2 ) - a(x 2 , (r - r') 2 ) . (2) 

After integration over Zj's and summation over n we obtain the following formula 

lUj " ApA) - C F x 1 G{x 1 ,ml) [d 2 r* G (h,r,z = 1/2) * v (r)® [ d 2 r'* G (h,r', z = 1/2) * v (r') 



dYd 2 b 

x i| exp ( - a (x 2 , (r - r') 2 ) p2iZ^ + exp ( - (<j(ai2,r) + <x(x 2 ,r / ) + a in (x 2} r,r')) p2R A S j 

-2cxp(-(<7(x 2 ,r)+a(x 2 ,r'))p2R A n . (3) 



The color factor in ([T]) and ([3]) is calculated in the N c ^> 1 limit. 

In the quasi-classical approximation the gluon saturation is given by 0] ^(a;) = 4 7r 2 a 2 pT(b), where p is the 
nucleon density in a nucleus, N c is the number of colours, b is the impact parameter and T(b) is the optical width 
of the nucleus. Q s .a determines the scale of the typical transverse momenta for inclusive gluon production [To| . Its 
value was extracted from the fit of the multiplicities of nuclear reactions at RHIC [ill, EH and from fits of the F 2 
structure function in DIS 0, Q3, QJ, [T3| . 

To generalize (|3|) to the case of AA collisions we use the approach suggested by Kovchegov [13] and replace the 
lowest order gluon field correlation function by the full MV formula as follows 
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Furthermore, noting that the dominant contribution to the integral over r' in Q comes from the region r > r' @ we 
obtain 



1 da (A A) 
S2 dYd 2 b 



'dCC 9 K 2 (C)e-^^ ixi) + Q ^ {x2) l (5) 



Eq. (0 is derived in the quasi-classical approximation which takes into account multiple scattering of the cc pair in 
the cold nuclear medium. At forward rapidities at RHIC and at LHC the gluon distribution functions ([4]) evolve 
according to the evolution equations of the color glass condensate. Inclusion of this evolution in the case of J/\V 
production presents a formidable technical challenge. Therefore we adopt a phenomenological approach of [l2| in 
which the quantum evolution is encoded in the energy/rapidity dependence of the saturation scales. 

In our numerical calculations we take explicit account of the impact parameter dependence of the saturation scales 
of each nucleus. We employ the Glauber approximation and assume that the nucleons are small compared to the 
size of the nuclei. The number of J/vP's inclusively produced in ion-ion collisions at given rapidity Y and centrality 
characterized by b reads 



dN AA (Y,b) 
dY 



C 



dNPP(Y) 
dY 



/ d 2 sT Al (s)T A2 (b-s) (Ql Al (x u s) + Ql A2 (x 2 ,b-s)) — 2 



f d( C 9 K 2 (C) exp (-^- {Q 2 SiAl ( Xl ,s) + Q 2 s , A2 (x 2 ,b- s)) 
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where X\/ 2 — m ''^' t e Ty with 
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-p\ , pt being the transverse momentum J/\V. The overall normalization 

constant C includes the color and the geometric factors C 2 ?/ {Air 2 a. s Sp) where S p is interaction area in proton-proton 
collisions. C also includes the amplitude of quark-antiquark annihilation into J/\E' and a soft gluon in the case of pp 
collisions. This amplitude as well as the mechanism of Fig. [T]A have a significant theoretical uncertainty. Therefore, 
we decided to parameterize these contributions by an overall normalization constant in ([6]). 
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FIG. 2: (a) J/ip rapidity distribution in Au-Au collisions for different centrality cuts, (b) Nuclear modification factor for J/S& 
production in heavy-ion collisions for different rapidities. Experimental data from [l9l ]. 



The rapidity distribution of J/tp's in pp collisions, the factor dN pp /dY appearing in Eq. [HI is fitted to the experi- 
mental data given in [1 81 ] with a single gaussian. In figure [2ja) the results provided by Eq. ([6]) are then compared to 
the experimental data from PHENIX Collaboration [l9[ for Au-Au collisions at y/s = 200 GeV. The global normal- 
ization factor C is found from the overall fit. There are no other free parameters. The agreement of the theoretical 
results with experimental data is reasonable. It is evident that the effect of the gluon saturation on the J/ip rapidity 
distribution in nucleus- nucleus collisions is to make its width a decreasing function of centrality. The distribution in 
the most central bin is significantly more narrow than in the peripheral bin. 

It is important that we describe well the data in the semi-peripheral region. This ensures that our model gives 
a good description of the J/\V production in dAu collisions. We also note that an earlier approach Q in which the 
same model was employed (albeit with an oversimplified nuclear geometry) provided a reasonable description of the 
data. Still a more detailed investigation is required which takes into account the exact deuteron and gold nuclear 
distributions. This will allow a model-independent fixing of the overall normalization constant C . Such an analysis 
will be presented in the near future. 
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Previously, the initial-state effects on the nuclear suppression of J/ip's have been estimated [26|, [2?J through the 
product of nuclear modification factors measured in dAu collisions. This estimate holds as long as the collinear 
factorization holds for the process of J/Psi production (this would be true for example in the limit of small cc-bar 
dipole size). Here we find that the effects beyond the collinear factorization are strong, and may account for a large 
part of the suppression measured in Au-Au, especially at forward rapidities (where the density of gluons in the initial 
state, inside one of the colliding nuclei, is the largest). 

To emphasize the nuclear dependence of the inclusive cross sections it is convenient to introduce the nuclear 
modification factor 

dN AA 

R A A(y,N pait )= dy dNPP . (7) 

iVcol] - 3 g- 

It is normalized in such a way that no nuclear effect would correspond to Raa = 1. In Fig. OR) wc plot the result 
of our calculation. The nuclear modification factor exhibits the following two important features: (i) unlike the open 
charm production, J/\& is suppressed even at y = 0; (ii) cold nuclear matter effects account for a significant part of 
the "anomalous" J/^P suppression in heavy-ion collisions both at y = and y = 2. 

III. CONCLUSIONS 

The main results of this paper are exhibited in Fig. [2j It is seen that the rapidity and centrality dependence of 
J production are reproduced with a reasonable accuracy even without taking into account any hot nuclear medium 
effects. This observation allows to conclude that a fair amount (and perhaps most) of the J/^ suppression in high 
energy heavy-ion collisions arises from the cold nuclear matter effects. In other words, J is expected to be strongly 
suppressed even if there were no hot nuclear matter effect produced. In this sense we can talk about the separation 
of the cold and hot nuclear medium effects as advertised in the Introduction. 

A certain fraction of J/vt suppression in forward direction can be attributed to suppression of the constituent c and 
c quarks [2CJ, |2l|, |22j, |23[ which would lead to suppression of D- mesons |24j. However, D- mesons are not predicted to 
be suppressed at central rapidities [24| ■ 

The reason for J/vff suppression at mid-rapidities is that the multiple scattering of cc in the cold nuclear medium 
increases the relative momentum between the quark and antiquark, which makes the bound state formation less 
probable. It was proven in [2o| that unless quantum log(l/a;) corrections become important, the inclusive gluon 
production satisfies the sum rule that requires the nuclear modification factor to be of order unity. Similar sum rule 
holds for heavy quark production but fails in the case of a bound states, such as J/^. 

We realize that although our calculation gives the parametrically leading result at high gluon density, other pro- 
duction channels involving the gluon radiation in the final state and the color octet mechanism of J production 
may give phenomenologically significant contributions. These are likely to become the leading mechanisms in the 
peripheral collisions where the strength of the gluon fields is significantly diminished. However, we believe that our 
main results are robust for central high-energy collisions of heavy ions. These results imply that the observed J/^J 
suppression is a result of an interplay between the cold and hot nuclear matter effects. 
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